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We study the non-equilibrium dynamics of an impurity in an hiarmonic trap thiat is kicked with a well-defined 
quasi-momentum, and interacts with a bath of free fermions or interacting bosons in a ID lattice configuration. 
Using numerical and analytical techniques we investigate the full dynamics beyond linear response, which 
allows us to quantitatively characterise states of the impurity in the bath for different parameter regimes. These 
vary from a tightly bound molecular state in a strongly interacting limit to a polaron (dressed impurity) and a 
free particle for weak interactions, with composite behaviour in the intermediate regime. These dynamics and 
different parameter regimes should be readily realizable in systems of cold atoms in optical lattices. 
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Impurities play a crucial role in determining the low-ternperature features of a number of condensed matter systems. These 
impurities may be localized ones, as for the x-ray edge and Kondo effects or mobile ones, like the itinerant single 
electrons modified by the phonon bath of the solid state crystal in which they move (polarons) IH, or the single spin-flipped 
electron moving in a lattice populated by opposite spin electrons, as studied in the context of high Tc superconductors 1^. With 
recent experimental developments for systems of ultracold atoms, such as the the tunability of the two-body interaction with 
the aid of Feshbach resonances 101, a particularly well controlled environment to explore the properties of these types of 
many-body systems has become available. 

For two- and three-dimensional systems, these advances have enabled the experimental ll6l- [Tlll and theoretical lfl2l - l23ll study 
of mobile impurities inside a fermionic bath, i.e. a type of polaron, usually created by preparing two-component ultracold 
Fermi gases with a large number imbalance between the components. In that context, the possibility of tuning the bath-impurity 
interaction across a wide range and even from attractive to repulsive regimes has opened up both the polaronic and the molecular 
regime to investigation. These advances have also stirred interest in using these fermionic systems to study the dynamics of the 
x-ray edge effect, which is induced by a localized impurity 112311 . 

At the same time, the ability to restrict the spatial dimension of the ultracold gas experiments almost arbitrarily has also made 
the study of impurities inside one-dimensional many-body baths possible. As movement of the impurity in such a bath can 
very easily involve the collective motion of many bath atoms, the result can be profoundly modified compared with what would 
be expected in higher-dimensional systems, giving rise to a regime of subdiffusive impurity motion, in which it can displace 
only proportional to the logarithm of time, slower than any power law ll24l - l28ll . Another reason for the particular interest in ID 
impurity-bath systems is that they make particularly compelling benchmark systems for a wide range of impurity-bath systems, 
due to the powerful theoretical approaches available to treat interacting ID svstem s ll29l[30ll . For example, the ground state of an 
impurity in a ID Fermi gas in a lattice was calculated via exact numerical methods l3l[l . demonstrating that it can be described by 
a polaron-type ansatz for weak interactions, while the strong interaction regime corresponds well to the strongly interacting limit 
of the Bethe ansatz. Static properties of polarons in ID ultracold Fermi gases have been studied also in ll32ll . and recently there 
has been interest in exploring the dynamics as well ll33l[34ll . Complementary to the fermionic case, the dynamics of an impurity 
in a continuous bosonic bath was studied recently experimentally and theoretically Issl - lssll . Major advances with single-site 
addressing and manipulation in optical lattices Il39ll40ll have recently enabled the reaUzation of lattice impurities within a bath 
described by a ID Bose-Hubbard model lliUl . 

In this article we explore the basic dynamical properties of a single impurity in a lattice potential and a harmonic trap in ID, 
which interacts with a bath of free fermions or interacting bosons, also confined in the ID lattice. Specifically, we consider 
the non-equilibrium response of the impurity to a kick with well defined momentum. A key open question in this context is 
how to characterise the role that the bath atoms play in the dynamics. In particular, we ask whether the dynamical response of 
the system implies polaronic behaviour, in which the properties of the particle are renormalised by the presence of the bath, or 
whether the interaction gives rise to other states, e.g., to tightly bounds pair or more complex objects. 

Using time-evolving block decimation (TEBD) methods ll42l - l45ll in conjunction with Bethe-ansatz results, we study the non- 
equilibrium dynamics of the impurity beyond the weak coupling assumptions of linear response theory. We show that the 
observed oscillation frequencies of the impurity-bath system can be mapped onto different physical states, and explain their 
dependency on bath density and strength of the impurity-bath interaction. In different limits we see that the behaviour ranges 
from a tightly bound pair for strong interactions to polaron-like behaviour at weak interactions for a fermionic bath. The latter 
case is characterised by an interesting internal dynamics corresponding to the scattering between a bound pair and an impurity 
particle propagating through the fermionic bath. We also compare these results with the case of a bosonic bath to better define the 
role of the Fermi sea. Generally, we find that the physics for a boson bath can be qualitatively and even quantitatively similar to 
the fermionic case, for both the doublon and the polaron regime, provided the boson-boson repulsion is larger than the attractive 
interaction between bath and impurity. 

This setup and characterisation of the dynamics should be readily realisable with cold atoms in optical lattices, and we expect 
our zero-temperature results to hold also at finite temperature, provided it is lower than the energy scale given by the oscillation 
frequency. 

This article is organised as follows: we first introduce the system and describe the method used in Section U There we also 
discuss the effects of combining a lattice potential with a trap in the case of an impurity that does not interact with the bath. 
In Section we present numerical results for the impurity dynamics and explain them both in the regime of strong and weak 
interactions through comparison with analytical methods. We especially discuss the frequency spectrum of oscillations, and 
identify different physical regimes of impurity behaviour In Section |IIll we compare with the case of a bosonic reservoir, and 
identify similarities and differences to the case of the fermionic bath dependent on the boson-boson repulsion. In Section |IV] 
we discuss our findings and make a connection to earlier impurity and polaron studies. Finally, an appendix contains details of 
several analytical results we have derived. 
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FIG. 1: Time dependence of the local doublon density {rii^riii) for \U\/J — 10 and N^= 20 (a), 60 (b), 140 (c), 180 (d). Horizontal axis: 
site index i, vertical axis: time (in units of J^^). In the simulations we first calculate the ground state of the lattice (number of lattice sites 
L = 200) loaded with one impurity (N_i = 1), in presence of a bath of spin-up particles (different simulations were performed with varying 
numbers of up particles Nf = [1, 10, 20 . . . 180, 190, 200]), whose description is provided by the Hubbard Hamiltonian. In addition to 
the attractive interaction between spin-up and spin-down particles, characterized by the parameter U, the impurity is confined by a parabolic 
potential of the form Vi^, with i the site index, and V/J = 0.1. Att = 0^, the down particle is kicked with quasi-momentum k — O.Itt. 
As is described in the text, it is possible to see that, as a function of A^^, the dynamics exhibits different regimes, characterized by a transition 
from doublon-oscillations to a regime where the internal dynamics of the polaron dominates, and eventually to free-particle oscillations. All 
simulations are performed with a TEBD code, with the following numerical parameters: lattice size L = 200, unless otherwise stated; Schmidt 
number x = 80; initial imaginary timestep (ground-state calculation) Sti — 0.1J~^; timestep (real time evolution) St — 0.02J~^. After 
testing with larger values of the Schmidt number, we have used the value x = 80 which provides accurate results both for the ground state and 
the time evolution. The reason of the effectiveness of this rather small value of x fies in the reduced size of the Hilbert space for a strongly 
imbalanced gas. 



I. FERMIONIC SYSTEM 
A. Basic model and method 



We consider a setup that is constituted by a optical lattice, loaded with a number-imbalanced mixture of two hyperfine species 
of fermionic atoms, hereafter labelled t and I, which are confined to move along one dimension. Our interest lies in the case of 
extreme imbalance, namely = 1 and iV-|- e {1 . . . L}, where L is the lattice size and N^/i is the total number of atoms for 
each species. In addition to the optical lattice, the I atomic impurity experiences a parabolic confining potential. For atoms in 
the lowest Bloch band, the system can be described by the Hubbard Hamiltonian 

H = -jy^c|^Q+ig + h.c. + U y^^nj^mi + Vy^^nii (i - ^ ^ j , (1) 

icr i i ^ ' 

where J represents the hopping amplitude between neighbouring sites, [/ < is the on-site (attractive) interaction energy and 
V characterises the strength of the parabolic confining potential for the impurity. Throughout the paper we set H = 1, and we 
choose as the length scale the lattice period a. Therefore, all energies are given in frequency units and all momenta are given in 
units of 1 /a. 

To obtain the ground state of this Hamiltonian and simulate the full many-body dynamics after the impurity has received a 
kick with a defined quasi-momentum, we use a code based on the TEBD algorithm, for which more details can be found e.g. 
in Refs. ll46l - l49ll . In our simulations we have considered a lattice size of L = 200 sites (L = 400 in one case), ~ 1, 
= [1, 10, 20... 180, 190, 200], and \U\/J = 0,0.5, 1, 3, 5, 10, 20 with particular emphasis on the case \U\/J = 10. At 
t = 0+, a quasi-momentum k = O.Itt is imparted to the l impurity. 
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FIG. 2: Time dependence of the local doublon density (nifnn) for \U\/J = 1 and N-^= 20 (a), 60 (b), 140 (c), 180 (d). For \U\/J = 1, the 
time evolution of {ni^riii) is characterized by single-particle oscillations decreasingly damped for increasing bath particle population. For 
A'^ = 20, small-amplitude oscillations (not visible here, see Fig. ^ appear; these oscillations are associated to pair-breaking mechanism as 
we will discuss in Sec. Ill CI All simulation parameters (except U) are the same as in Fig. [T] 



B. The non-interacting impurity 



To provide a better description later on of the effects on the motion of the impurity that are induced by interactions with the 
bath, we first single out the dynamical effects induced by the concomitant presence of the harmonic trapping and the ID lattice 
aiU — Q. That is, we consider a single particle in a potential formed by combining a lattice and a harmonic potential. As shown 
in ifsoll . if the particle is initially in the ground state of the combined potential, and then the harmonic trap is displaced by an 
amount 5, the expectation value of the particle position is 



{x) = 6 exp 



sin^ [V/U) 



cos 



V 
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Ho 



sm 



Vt 
T 



(2) 



where V is the strength of the trapping potential, and moff is the effective mass of the particle in the lattice, which for our choice 
of units { a = 1) is = 0.5J~^. The harmonic oscillator length is defined via the trapping frequency lu and the mass of 
the particle as Oho = l/V^cffW, and on the other hand the strength of the harmonic oscillator potential V — mcguj'^/2: thus 
flho = (2TOoffy)~^/^ and with moff = 0.5J~^ we then should have a^o = {V/ J)^^^^. 

Our initial conditions involving a finite quasi-momentum kick differ from those in Ref. ifsoll , where the harmonic trap is 
initially displaced by an amount 6. To account for this difference, we need to introduce a 7r/2 phase in the cosine term of Eq.|2] 
and also compute the value of 6 corresponding to our initial momentum kick. We can do this by matching the energy of the 
initial state in each case. The e nerg y for the state at t — 0+ is given by the lowest eigenenergy of the combined harmonic trap 
and lattice system (eq. (14) of ||5(||), with the addition of the kinetic energy given by the kick of a quasi-momentum k. If, by a 
semiclassical argument, we assume that this energy must be converted to potential energy (removing the zero-point energy from 
both terms), then we have 
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(3) 



The left hand side of Eq. ([3]) represents the small V expansion of the potential energy of a particle in a ID lattice in the presence 
of harmonic confinement when displaced by 6 from the minimum of the potential, while the right hand side is its kinetic energy. 
We can thus deduce an approximate expression for S 
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Vmes 16 
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(4) 
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Combining Eqs. and (01, we obtain the relation between the centre of mass oscillation frequency ujcom and the theoretical 
value of the mass of the different particles 



As a benchmark, we have compared the frequency given by Eq. (|5]l with the value obtained from the numerical simulations for 
a free particle (nics — 0.5), and we have found good agreement between analytical and numerical values. 



In the following, we characterize the different physical states which the impurity may form inside the bath. The key ob- 
servable in our analysis is the oscillation of the time-dependent doublon density, defined as {ni\ni{){t), after a kick with 
quasi-momentum k = O.Itt has been imparted to the impurity at t = 0+. We find that this is a more useful quantity than 
the impurity density {ni{){t), as it provides more direct information corresponding to the bath dynamics and pairing with the 
impurity. 

Examples of the oscillatory motion of the doublon density are shown in Figs.[T]and|2]in the regime of strong {\U\/ J ~ 10) and 
weak (\U\/ J = 1) attraction respectively. In the strongly attractive case, we see the oscillation frequency shifting as a function 
of the bath filling, increasing only slightly while the bath density is below half filling, but jumping to much larger values above, 
from where it decreases again as the density approaches integer filling. Conversely, in the case of weak attraction the oscillation 
frequency only decreases slightly for low and high densities, and is roughly the same for all intermediate densities. 

The observable that encapsulates all these behaviours is the doublon centre of mass {X^]), defined as 



which is extracted from the full density and is shown in Figs.[3]and|4] 

As is explained in detail in the following sections for the different interaction and bath density regimes, we gain insight into 
the physics of the system by analysing the Fourier transform of X^i_{t), A'^|(a;), which is shown in Figs. |5]and|6]for strong 
and weak attraction. For strong attraction, X^^-iluj) shows that the oscillation of a tightly bound on-site pair dominates the 
dynamics at low density, while a polaron-like state is present but weak (low- and high-frequency peaks in Fig. |5}a respectively). 
The relative weights of the polaronic and bound-pair peaks reverse as the density of bath atoms increases above 0.5, with the 
polaron component becoming predominant and the bound-pair peak almost vanishing (high- and low-frequency peaks in Fig. 
|5}b respectively; notice e.g. how for ~ 140 the bound-pair peak has become essentially just a broad shoulder). 

In order to understand the physics of the bound pair for strong interactions between impurity and bath atoms, we make use 
of the so-called string hypothesis from the Bethe-ansatz solution of the Hubbard model (c.f. Appendix IV Al ). Using this, we 
see that a tightly-bound on-site pair dominates at low bath fillings, and can be understood in a two-body picture, where we can 
compute formulas for the pair oscillation frequency, as shown in Fig. [T] On the other hand, the polaron-like state dominating at 
higher densities can be explained by an essentially free particle oscillating inside a cloud of bath particles of size C^. This cloud 
forms due to the impurity acting very similar to a hard-wall boundary condition on the bath particles at its position - the density 
of bath particles around the impurity behaves like Friedel oscillations in a gas of free fermions at a hard wall. Subsection III Bl 
describes this in detail, including the change in the polaron oscillation frequency with density. 

By contrast, in the case of weak interaction, a polaron component to the oscillation is significant only for densities at or 
below 0.2 (low-frequency peak in Fig. |6}a), while the dominant component of X^i{ijj) stems from the motion of a free particle 
(low- and high-frequency peaks in Figs. |6}a and b, respectively). In this regime the polaronic component corresponds to the 
resonant scattering between a bound pair and an impurity particle propagating through the background particle bath. This will 
be described in section ilTCl 



In a ID Hubbard model with strong interactions, Bethe-ansatz techniques provide two sets of solutions, which together cover 
all possible eigenstates of the system - provided we assume the so-called string hypothesis to be correct Isoll . In each of these 
two sets of solutions, A and B, the total energy E and quasi-momentum P of a many-body state can be expressed in terms of 
Bethe-ansatz quantum numbers {/s;} and A. 




(5) 



II. DYNAMICS OF AN IMPURITY IN A FERMIONIC BATH 



A. Strong interactions: the low-frequency pealc 
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FIG. 3: The doublon centre of mass ^ii) for j(/|/J = 10 and (a) N-^ = 20 (purple curve), 60 (green dash-dotted curve), 100 (black 
continuous curve); and (b) A'^^^ = 100 (black continuous curve), 140 (blue dashed curve), 160 (red curve). In panel (a), it is possible to see 
how the centre-of-mass oscillations associated to the doublon dynamics are increasingly damped when the bath population increases towards 
half filling. Above half filling, the oscillations associated to the polaron internal dynamics (described in Sec. IIIBb start to appear, reaching the 
free-particle oscillation frequency for 7V-|- — >■ L. The deviation of X^^ for t — > lOOJ^'^ is a finite-size effect. It corresponds to the time for 
which excitations in the bath are reflected back from the system boundaries to the centre of the trap, affecting thus the doublon dynamics. All 
simulation parameters are the same as Fig.[T] 



For the solutions of type A 



N N 

P = ^ kj E = — cos [kj) + const. 

3=1 J=l 



(6) 



and for the solutions of type B, 



N-l N-1 

P = ^ kj + 2q E = — 2J ^ cos {kj) — 4J cosh ^ cos (g) + const. 
i=i 3 



(7) 



where ^ is defined by 



cosh ^ = W 1 



C/2 



16 J2 cos^ [q) 



In the strong coupling limit, it is possible to show that 



4Jcosh^cos(j -^U - 4JV \U\ - 4JV \U\ cos {k) + O 



C/2 



(8) 



(9) 



where k = 2(7 is the total quasi-momentum of the pair 

The eigenstates of type A correspond to the effectively free motion of both bath particles and the impurity, while solutions of 
type B correspond to a bound state of the impurity with one bath atom (with the remaining 7V^— 1 bath atoms again moving freely) 
(see Appendix IV At . For an attractive interaction, the B-type solutions are always energetically favorable, and for \U\/ J > 2 
there is no overlap between the bands associated with type-A and type-B solutions. 

Thus focussing on the B-type manifold of eigenstates, we can describe the observed behaviour of the low-frequency peak in 
X^l{uj) at large \U\ (c.f. Fig. |5]l by deriving the explicit expression for the many-body energy E in Eq. (|2l) as a function of the 
doublon quasi-momentum k (see Appendix lV Aj to be 



E^ = ^^yU^ + lGJ"^ cos2(k/2), 



(10) 
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FIG. 4: The centre of mass Xi^i{t) for |f/|/ J = 1 and (a) A^^ = 20 (purple dashed curve), 60 (green dash-dotted curve), 100 (black continuous 
curve), and (b) A''^ — 100 (black continuous curve), 140 (blue dashed curve), 160 (red curve). For \U\/ J = 1, the oscillations are dominated 
by the single-particle oscillation frequency co/ J ~ 0.62, and increasingly damped for decreasing bath population. All simulation parameters 
are the same as in Fig. |2] 




FIG. 5: The doublon centre of mass Af-f 4_(aj), \U\/J — 10. The Fourier analysis allows a quantitative estimation of the oscillation frequencies 
for different values of Nf, complementing the description given in Fig. [3] (a) A^^t = 20, 60, 100 (purple dashed, green dash-dotted and black 
continuous curves). For Nf = 20 only the doublon peak for oj/ J ~ 0.5 is present, while for Nf = 100 the polaron internal dynamics peak 
appears, for uj/J ~ 2.1 respectively, and for A^^ = 60 the precursor of the polaron peak starts to appear around u/ J ~ 4.1 . (b) A^^ — 
100, 140, 180 (black continuous, blue dashed, and red dash-dotted curve). The polaron internal dynamics peak shifts to the left for increasing 
interaction {uj/J ~ 2.1, 1.1, 0.8 for N-f = 100, 140, 180), while the doublon peak is increasingly damped. All simulation parameters are the 
same as in Fig. [T] 
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FIG. 6: X^iiuj), \U\/J = 1. (a) = 20, 60, 100 (purple dashed, green dash-dotted and black continuous curves) and (b) iV^ = 100, 140, 
180 (black continuous, blue dashed, and red dash-dotted curve). For Nf =20 both the polaron and the free particle peaks are present (oj/ J ~ 
0.4, 0.8 respectively, while for larger values of A'^, only the free particle peak is present. All simulation parameters are the same as Fig. |2] 



where the dependency of the energy on the other iV^ — 1 quasi-momentum quantum numbers { ki } has not been taken into account, 
being irrelevant in the dynamics here considered. From this, the effective mass of the tightly bound pair can be extracted to be 



^doublon 



1 



16.P 



c/|/,/>i 



4J2 



(11) 



which, in the limit \U\/J ^ 1, conincides with the expression that would be obtained from second-order perturbation theory. 

It is interesting to note that our Bethe ansatz solution delivers the same result for the quasi-momentum dependence of the 
doublon contribution to the energy as the simple solution to the problem of two distinguishable bound particles on a lattice 
does lf5lll52ll . In Fig.|7]we show the comparison between the numerical value of the doublon peak oscillation frequency (low- 
frequency peak in the strong-interaction regime), and the theoretical value obtained using our value for the doublon mass (fTTT i in 
conjunction with Eq. (|5]i. The excellent agreement shows that the doublon dynamics is well-captured by our Bethe-ansatz based 
model. 



B. Strong interactions: the high-frequency peak 



In the case of strong attraction (large \ U\) treated in this section, the oscillations of the doublon density also show another, high- 
frequency component in X^i{t), which is weak for low bath densities, but which comes to dominate the oscillatory behaviour 
above half-filling (see Fig. |5]l. Here, we argue that this feature can be understood in terms of the oscillation of the impurity inside 
the cloud of bath particles that surrounds it and which may restrict its motion when the initial displacement of the impurity by 
the kick at f = 0+ is larger than the size of this cloud. 

To describe the behaviour of the high-frequency peak in the following in more detail, we first characterise the size of the 
modified area for the static impurity, and then consider its effect on the kicked impurity. Since the polaron is a concept associated 
with the impurity surrounded by a distorted cloud, we will, in a loose sense, call the size of the density modified area the "polaron 
size". 



1. Effective size of the bath area modified by the impurity 



With attractive interactions and an impurity that is localized by the tight parabolic confinement potential, the density of bath 
particles is modified locally in the center of the system. Let us denote the size of this modified region of the bath by C,. We make 
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FIG. 7; Comparison between the numerical value of the low-frequency peak oscillation frequency for the large interaction case (as extracted 
from X^i{u)), and the theoretical value for |I7|/ J > 1, given by Eqs. JlOt Jilt). This comparison allows us to establish that the low-frequency 
oscillations of 4, can be described in terms of oscillations of a particle (composed by one spin-up and one spin-down particle) in a parabolic 
confining potential. All simulation parameters are the same as in Fig. [T] except the interaction strength and the lattice size. Here \U\/ J =3, 
5, 10, 20 and L = 200 or 400. The error bars account for the finite resolution of the Fourier transform. Larger lattices allow for longer 
simulation times before the reflection from the edges starts playing any role in the doublon dynamics, thus allowing higher resolutions in the 
frequency domain. It is worth noting that the full Bethe-ansatz formula works even for moderately weak interaction. The characteristic size of 
the composite particle is given by (see Fig.|15ll, corresponding to a value of ~ 2.3 for \U\/ J — 3. 



theoretical solution 




a first estimate for the value of ^ by considering the limiting case U —00, with the impurity being located on a single site, i. 
In this case, due to the limit of infinite attraction, the site i will act as a hard- wall boundary condition llssilsill . and the density 
profile of the bath particles around the impurity will undergo Friedel oscillations ifssll (see Fig. [8]l, as was recently pointed out 
in ll28ll for repulsive interaction and in absence of the lattice, according to the following formula 

We note that in a ID system the period of these oscillations is always independent of the strength of the interactions between the 
bath particles and the bath impurity interactions, and is always 2kpf for a fermionic bath, and 2ttpq for a bosonic one - all that 
changes with the interaction is the amplitude of the Friedel oscillations. 

In Fig. |9] we show how approximating either the on-site bound pair (in case of strong attraction) or the single particle (for 
weak attraction) as a boundary condition located at the minimum of the parabolic potential works well to describe the density 
modulation induced in the bath when is large, while it fails in the case of a more dilute system. As we will use this 
approximation only in the regime of intermediate and high filling, the discrepancy at low filling is not a problem. 

In the following subsection, we will assume that these Friedel oscillations surrounding the impurity can be approximated by 
a box-like area, inside of which the impurity can move effectively like a free particle, regardless of interaction strength. 

We assume that the effective box size is given by C = aCm, where is the position of the first extremum of Pf{j), and a < 1 
is a constant determined by fitting to the numerical data. The position of the first extremum of pf{j) is approximately given by 
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FIG. 8: Friedel oscillations (see Eq. J12b ) induced by the presence of a hard-wall boundary condition at j = {L — l)/2 for (a) 20 (purple 
dashed curve), 60 (green dash-dotted curve), 100 (black continuous curve), and (b) 100 (black continuous curve), 140 (blue dashed curve), 
180 (red dash-dotted curve) bath particles respectively (see Eq. l ll2t ). 
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FIG. 9: Comparison between the numerical estimate of the ground-state spatial oscillations induced in the bath by (a,c,d) a single down 
impurity or (b) a doublon (blue continuous curves), with the Friedel oscillations induced by hard-wall boundary conditions (green dash-dotted 
curves), with (a) iV^ = 180, \U\/J = 10, (b) = 20, \U\/J = 10, (c) = 180, \U\/J = 1, (d) iV^ = 20, \U\/J = 1. In the 
high-density limit, the agreement between the numerical calculations and the approximate model is good, while, in the low-density limit, the 
numerical results suggest that the role of the impurity cannot be described in terms of an impurity-induced Hartree potential. This hypothesis 
is confirmed by the exact diagonalisation of an Hamiltonian describing a spin-polarized gas in presence of a localized potential mimicking the 
potential induced by the impurity. In this case the results are compatible with the results implied by Eq. l ll2t . All simulation parameters are 
the same as in Fig. [T]and Fig. |2] 



2. Description of the dynamics in terms of the effective modified area 

Using the concept of effective modified area (polaron size) introduced above, we can now explain those impurity dynamics 
which lead to the high-frequency component in X^]^{t). There are two qualitatively different regimes. In one regime, C is smaller 
that the displacement 6 (c.f. Eq. (|4]l) which the atom experiences in the combined trap and lattice potential due to the initial kick. 
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FIG. 10: Schematic depiction of the effective modified area effect (polaron size), describing the high-frequency peak for the strong-interaction 
limit. The doublon (black arrows) is restricted to oscillate within the area defined by (red arrows). Depending on the magnitude of the initial 
displacement of the impurity, S, the impurity can exhibit two different regimes of motion: For < S, the oscillation frequency of the polaron 
mode must adapt to the effective polaron size, leading to larger frequencies as <^ shrinks. Once the density drops sufficiently for ^ > 5 to hold, 
the frequency of the polaron-mode shifts closer and closer to the frequency of the free particle. See also Fig. [TT] 



In the other regime, ^ is larger than 5. 

For C < the dynamics of the impurity is constrained within the modified area of size ~' (, c.f. Fig.[TO] In a non-interacting 
system, the kick would cause the impurity particle to be displaced by the full amount of S, but in the presently considered 
regime of large \U\, the cloud size C, actually constrains the motion. This confinement can be understood in terms of energy 
conservation: it is energetically unfavourable for the impurity to exit the cloud of bath particles it has attracted around itself in 
the initial ground state, a mechanism analogous to that of repulsively bound pairs ifsill . The confinement, however, is not perfect: 
part of the energy and quasi-momentum leaks out from the boundary given by ( during each oscillation cycle, thus leading to 
damping. 

The oscillation frequency of the impurity in the box can be obtained by a semiclassical argument, considering a particle 
bouncing off rigid walls placed at a distance 2(, leading to the following expression for a;,; 

„ fc/m k L- ^ 

c^C = 27r4— =^ -^J, (15) 

% a L 



for a quasi-momentum kick k, and a particle of mass m IdSII . Comparing Eq. ( fTSl ) to the TEBD-computed frequencies of the 
high-frequency modes of X-^^{t), it is possible to deduce a value for a ~ 0.25. 

On the other hand, for C > 5 the oscillation frequency is determined by the natural oscillation frequency of a free particle in 
an empty lattice with a superimposed parabolic potential: the constraint does not play any role since the motion of the (free) 
particle takes place in a region which is smaller than the size of the induced bath cloud. 

Fig, nn summarises our findings, showing that with a fitting factor a ~ 0.25 for the effective size of the bath density cloud, 
our model is able to explain the shift of the high-frequency mode in the center-of-mass motion of the doublon Xf i with density. 



C. Weak interactions 



As anticipated, for interactions \U\ < 1, the peak in X^i{oj) associated with the oscillation of a tightly bound on-site pair 
is not present. Nevertheless, we observe two distinct modes in the doublon-density center-of-mass oscillations X^i{t) in the 
weakly interacting regime as well. One of them appears due to the free oscillatory motion of a non-interacting particle as derived 
in Sec. IIBI Another peak, at lower frequencies, appears as well (see Figs. l6land[T2]). 

The underlying physics of this low-frequency peak derives from a resonant transition between a bound pair and scattering 
states, specifically the spin-down impurity at zero quasi-momentum and a spin-up bath particle at the Fermi quasi-momentum 
kp^. The frequency of this peak can thus be obtained considering the difference between the energy of the (weakly-bound) pair 
and the energy of the scattering state 

ujpoi = + E^- E^i = -2 J [1 + cosikpt)] + Vu^TTeJ^. (16) 

Figure [12] shows the excellent agreement between this model and the numerical results from the TEBD calculations, within the 
error bars set by the finite resolution of the Fourier transform. 
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iV,./L 

FIG. 11: (a) Oscillation frequency (for the high-frequency peak in the \U\/ J = 10 case, corresponding to single particle rather than doublon) 
obtained from the simulations, ujnum, is compared to the non-interacting particle frequency, ojho, and to the frequency given by the effective 
confinement to the area of size C,, ujq. The plot shows the crossover from the latter to the former frequency for increasing bath population 
N-\ /L. (b) The crossover happens when starts to be considerably larger than the kick-induced displacement 5 in the non-interacting system. 
All simulation parameters are the same as Fig. [T] The errorbars account for the finite resolution of the Fourier transform. 
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FIG. 12: Comparison between the numerical value of the low-frequency peak at |[/|/ J = 1 and ujpoi corresponding to the resonant frequency 
of the bound-doublon to scattering-states transition. This peak can be seen in the A*'^ = 20 plot (small peak around uj ~ 0.2 in Fig.[6j. Above 
Nf/L ~ 0.2 it merges with the free-particle peak, and eventually disappears. In this picture, the spin-up bath particles are resonantly bound 
to the spin-down impurity, analogous to the case of static polarons. All simulations parameters are like in Fig. [2] The errorbars account for the 
finite resolution of the Fourier transform. 



Further, the existence of the two peaks in ^"-1-4,(0;) can be related to the known structure of the polaron ground state: in the the 
polaron ansatz |fl2| of the type 

\^)^VZcl\FSm^+ 'l^'^At^,tcUi\FSh\0)i (17) 

k>kl,,q<kl. 

the first term describes the impurity at rest in the presence of the unperturbed Fermi sea of the bath (with the quasiparticle 
weight Z), and the second term a coherent superposition of states, in which the motion of the impurity is correlated with a 
single particle-hole excitation of the Fermi sea (like the weakly bound state entering Eq. (fT6b). When such a state experiences 
an initial quasi-momentum kick att = 0+, both parts of Eq. iT7\ will contribute to the dynamics. Based on our results, we 
can hypothesize that the first term of Eq. (fTTI i leads to the high-frequency "free particle" peak in ^"^4^(0;), while the second 
term causes the low-frequency peak, which is caused by breaking up the correlated state between impurity and a particle-hole 
excitation of the Fermi-sea, as we have shown through the agreement of Eq. ( fT6] l with our simulation data (c.f. Fig. \T2\ . The 
relative height of the peaks could perhaps provide a useful tool for determining the relative weights of the two terms in Eq. ( fTTI i. 
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D. Damping 



Along the lines of the above discussion about the effective area induced by Friedel oscillations, it is possible to intuitively 
understand the oscillation damping. When approaching half-filling, both from the low- and the high-density limit, the oscillation 
damping is increased. This increase is associated to the decrease of the effective area: as we have shown above, the dou- 
blon/particle (low- and high-density limit respectively) are free to oscillate in the region |a;| < Close to half-filling, C 
and the energy initially transferred to the impurity is, in turn, transferred to the bath particles, since the (coherent) evolution of 
the doublon is confined to a region whose size is smaller than allowed by energy conservation. In other words, excitation of the 
oscillation does not conserve energy, thus the oscillation can only appear as a transient and the energy given by the kick has to 
be dissipated through the bath particles. Quantitative analysis of the damping is a subject of further study. 



in. THE KICKED IMPURITY IN A BOSONIC RESERVOIR 



Comparing the results obtained for a fermionic reservoir to those from a bosonic one enables us to state which features of 
the observed dynamics are universal, and which are particular to the fermionic bath. Towards this, we have performed TEBD 
simulations for a two-species Bose-Hubbard Hamiltonian 

=-J^&L&^+ia + /i.c. + C/^?i,tn,; + Y^n,t(n»t-1) + ^5Z"U («-^^-^) , (18) 

i a i i i ^ ^ 

where trap parameter V, tunneling J and the values of U and A'^ are the same as described in section U as is the impurity 
preparation and the initial kick to the impurity (as in the fermionic case we will set J = 1, a = 1). The two key differences are 
that now bier, are operators for softcore bosons, which interact repulsively on-site with energy W, if a =\. Here, we have 
considered W/J = 4, 10, 20. 



A. The weak interaction limit 



Repeating the TEBD simulations for weak interactions between the impurity and a bosonic bath, we find that the two-peak 
structure discussed in 111 CI persists for all values of W we have simulated, as shown in Fig. [T3] Moreover, the position of the 
polaron-dissociation peak (c.f. IllCb is still described remarkably well by the theory for the fermionic bath, even for the lowest 
value of W , W/ J = 4. These findings can be understood analytically by observing that, at low densities and reasonably large 
values of W, one-dimensional lattice bosons map to spinless fermions with weak nearest neighbor attractions. This mapping 
is achieved by describing the sector of low-energy, long-wavelength excitations of the ^ component of the Hamiltonian ( fISl l 
as a Tomonaga-Luttinger liquid (TLL), whose properties are characterized by the so-called TLL parameters, Ki, and vi, |5g]. 
Expanding the ^-sector in Hamiltonian ( fTSl l and the bath density operator n^f in terms of the canonically conjugate TLL field 
operators (f>-f{x) and Of{x), one obtains 

g2im[7rpoani-0t(an)] j ^ 



where expanding the impurity-bath coupling in this way presupposes that the impurity does not distort the bath density 
too much locally. In the limit of large W/J, the TLL parameters are known perturbatively, ~ 1 + simrpo, 

Ub = 2 J siuTrpo (l - ^Po cosTrpo) 

Now, a bath of spinless lattice fermions with attractive nearest-neighbor interactions Vnn and Fermi quasi-momentum kp 
coupled to an impurity with an on-site density-density interaction, can be mapped to a TLL in ID in a manner identical to ( fT9] l, 
where Kb is replaced by Kj and Ub by uj ll29ll . These parameters are also known from perturbation theory for small Vnn- 

Computing Kb and Ub for densities between and 0.2 and large W/ J shows that Eq. ( fT9l ) can be read equivalently as 
the model of an impurity coupled to weakly nearest-neighbour attractive spinless fermions. For example, for W/J = 10 and 
= 0.1, Kb = 1.039, Ub/J = 0.59, values that are best matched by Kf and Uf for |V^|/J = 0.75, kp = O.Itt. Crucially, 
a value of Kf ~ 1.039 is still very close to the values for free fermions, Kf ~ 1. Thus, the continued applicability of Eq. 
(fT6] l - which had initially been developed from a one-body picture of the free fermion bath - to predict the polaron peak in the 
weak-coupling regime even for a (sufficiently repulsive) bosonic bath, can be explained (c.f. Fig. [T3] ). 
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FIG. 13: Comparison of Xfi{Lj) between fermions and bosons in the weak-attraction case, |(7|/J ~ 1, A'^ ~ 20. The position of the low- 
frequency peak for the fermionic bath, which stems from the break-up of correlated states between the impurity and particle-hole excitations 
in the bath (c.f. IllCt . is very similar for the bosonic bath, at all values of W. See text for details. 



B. The strong interaction limit - higher-order bound states of the impurity. 

In the opposite limit of large bath-impurity attraction, \U\ ^ J, the kick-induced dynamics of the impurity start to depend 
crucially on the ratio As an example, as for the fermions, we focussed on the case \U\/J = 10. As long asW > \U\, 

the spectrum of doublon dynamics, Xfi{uj), remains qualitatively unaltered from the case of the fermionic bath: when < 0.5, 
the dynamics of the kicked impurity are dominated by the oscillation of the doublon-mode, whereas for > 0.5 this doublon 
peak increasingly flattens out and eventually disappears as n-f increases above half-filling, as shown in Fig. [14] At the same 
time, like for fermions, a high-frequency peak appears for > 0.5, increasing in amplitude as n-f grows above the threshold 
while the doublon peak decreases, signaling the transition of the dynamics to a regime where the unbound impurity oscillates 
inside an effective modified area (c.f. lIlBb . Even for a substantial value of W, W/ J = 20, the effective area is clearly smaller 
than in the fermionic case for a range of densities, as evidenced by an upshift in oscillation frequency over the fermionic case at 
equal density, and only becomes comparable to the size of the effective area for Fermions when nf w 0.9. 

On the other hand, when the boson-boson repulsion W becomes comparable to or smaller than the magnitude of the boson- 
impurity attraction \ U\, the numerics clearly show that higher-order bound states between impurity and bath particles are formed 
in the ground state. At W/ J = \U\/ J = 10, both doublon and trion states (the impurity binding to one or two bath particles 
on-site respectively) are present, whereas, for W/ J = 4 doublon, trion and quatrion bound states are occupied, with the trion 
state carrying the largest weight at any bath density. When quasi-momentum is applied to the impurity by the kick, these higher 
order bound states perform oscillations, at frequencies significantly lower than those for the doublons due to the even higher 
effective mass. Interestingly, the damping we observe becomes gradually smaller the smaller W/ J is, with the oscillations at 
W/ J = 4 showing almost no decay at any value of in the time domain over which we simulate. Whether this effect is due to 
the partial ability of ID superfluids to be protected against excitations Il29i, i57ll - which would be the source of any damping of 
the bound state oscillations - is an interesting question for further study. 
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FIG. 14: Comparison of Xfi{ij) between fermions and bosons in the strong-attraction case, \U\/J — 10. a) A'^^ — 140. For the fermion 
bath (black solid line), the peak from the unbound impurity oscillating inside the effective modified area (c.f. IIIBb has become dominant, with 
the doublon peak at lower frequency having almost completely flattened out. For bosons at W/ J = 20 (blue dashed line), the behaviour is 
qualitatively similar, but the oscillation of the doublon is still dominant, as evidenced by the larger low-frequency peak. For W/J = 10 (green 
dotted line) and W/J = 4 (red dash-dotted line), only the high-order bound states oscillate, there is no dynamics of an unbound impurity 
oscillating inside an effective modified area anymore, b) A'^ = 180. For the fermion bath, the peak from the unbound impurity oscillating 
inside the effective modified area (c.f. Ill Bt is now fully dominant, with the doublon peak at lower frequency having completely vanished. The 
dynamics for the boson bath at W/J = 20 (blue dashed line), now almost completely matches that of the fermions, showing that they also 
form the effective modified area like the fermions do. For W/ J = 10 and W/J = 4, only the high-order bound states oscillate as in a), there 
is again no dynamics of an unbound impurity oscillating inside an effective modified area. 



IV. DISCUSSION AND CONCLUSIONS 



The dynamics of the impurity moving on a ID lattice inside a fermionic bath, or a strongly repulsive bosonic one, shows 
intriguingly complex dynamics, as can be seen by studying the time-evolution of the doublon density. One of the characteristic 
frequencies that appears corresponds to the motion of a free particle, appearing in the limit of high filling of the bath and is 
easily understood due to the increasingly uniform interaction energy the impurity experiences on all sites. Another feature of 
the dynamics, namely the oscillations of a bound pair are also rather intuitive to understand, allowing to draw an analogy to the 
molecule vs. polaron question in three dimensional continuum systems. The bound pair is present in regime of large \ U\/J and 
low density, like the molecule in the polaron vs. molecule analogy. In the limit of small \U\/J,we observed the dynamics of a 
free particle side-by-side pair breaking of the correlated states of a polaron. This actually corresponds well to the polaron ansatz 
lfT2ll which is a superposition of a non-interacting Fermi sea (free particle) plus a contribution from correlated particle - hole 
states. Thus, the crossover from a polaron to a bound pair with increasing interaction is also taking place in analogy to higher 
dimensional continuum systems. In our system, however, we have a feature that does not have any analogy in the polaron vs. 
molecule crossover in 3D continuum, namely, the high-frequency peak in the large \ U\/J regime, which becomes dominant for 
bath fillings above 0.5, which is the result of the impurity oscillating inside the cloud of bath atoms dressing it. This process can 
be read as a kind of internal dynamics of the polaron. Observation of the dynamics predicted here should be feasible in currently 
available ultracold gases systems, provided that the temperature is below the energy scale of the oscillations which we found to 
be of the order of O.IJ-IJ. 



\U\ range 


Bath population 


Dynamics regime 


Strong interaction 


Large A^^ 


Free particle 


Strong interaction 


Intermediate Nf 


Bound pair -I- polaron internal dynamics 


Strong interaction 


Small iV-,. 


Bound pair 


Weak interaction 


Large A*"^ 


Free particle 


Weak interaction 


Intermediate & small A'^^ 


Free particle -I- polaron 



TABLE I: Table summarising the different regimes identified for the problem considered here. 
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V. APPENDIX 

A. Effective mass from the Bethe ansatz 

We show here how to gain some insight into the problem through the string hypothesis for the solution of the Lieb-Wu 
equations, whose solutions describe (most of) the eigenvalues of the Hubbard Hamiltonian in one dimension (Chapter 4 of 13^ ). 
in the limit of large lattice lengths L. For a fixed total number of particles N and number of down particles M, the patterns of 
which the solutions of the Lieb-Wu equations are composed can be classified in three different categories; 

• fc — A strings; 

• single real values of kj; 

• A strings. 

Every eigenstate of the Hubbard Hamiltonian can be represented in terms of a particular configuration of strings, containing Af„ 
A-strings, fc — A strings of length n (in our case n = 1), and Aie single kj. Here A/„, A/^, Aie are related to the N and M 
by 

OO 

M = J2n {M„ + A/,;) (20) 

n=l 

QO 

N = Me + ^'2nMl^. (21) 

ri=l 

In our case A/ = 1 implies the existence of two classes of string solutions 

A) Mn = 1, = 0, Aie = L: this solution is characterised by one A string constituted by a single real value, N real fc^s 
and no fc — A strings. 

B) A/„ = 0, A'jf/j = l,A^e = L ~ 2: in this case the solution is characterised by iV — 2 real fc^-s and one fc — A string, 
characterised by two (complex-valued) fci.2S and one real A, related by 

sin(fci,2) = A ± iU/iJ. (22) 
In terms of the string parameters, energy and quasi-momentum are given by 



P 



N-2M' OO M'r, 

^ kj - ^ ^ (2 Re arcsin (A'„" + niU/U) - (n + l)7r) 

j — l n— 1 a— 1 

N-2M' OO M' 



mod 27r (23) 



E = -2J cos(fcj) +4J^ ^Rcy 1 - (A'„ +™C//4J)^ (24) 

j—l n—1 a— 1 

For the two classes of solutions previously identified, see A and B above, take the following form 
A) 

N 



P = J2k, (25) 

N 

E = -2 J 'Y COS (kj) + const. (26) 
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FIG. 15: Pair size ^) as a function of U. 



B) 



N-l 



P ^ % - 2 Re [arcsin (A + iU/iJ)] 



E = -2 J ^ cos {kj) + 4JRc y 1 - (A + iU/ijf + const. 



(27) 
(28) 



The solution containing the fc — A string can be written in a more transparent form as 

N-l 



Af-l 



E = — 2J cos (fcj) — 4Jcosh^cos (q) + const. 



(29) 



(30) 



where g = Re [fci] = Re [fci] and ^ — Im [/ci] 
is possible to prove (see OOll p. 134) that 



Im [k2], with ki and k2 belonging to the fc — A string. For the fc — A string it 



cosh = A 1 



C/2 



16 J2 cos-^ (q) 



(31) 



From the form of the wavefunction associated to the fc — A string, ^ describes the (exponential) spatial decay of the pair, and 
thus can be interpreted as the size of the pair. 
In the strongly interacting limit we have 



4JcoshCcosg^ ;7-4JV|C/| -4J7|;7|cos(2g) + ( 



(32) 



which is coherent with the strong coupling calculation leading to the Heisenberg Hamiltonian for the Hubbard Hamiltonian in 
the strong coupling limit (modulo a [/ — > — [/ mapping). 

The spectrum is thus constituted by a lower band (fc — A+ single fcs solutions) and a higher band (single fcs and single A 
solutions). Intuitively, the former correspond to a pair in an unpaired background Fermi sea, while the latter to scattering states. 
The down particle which is kicked in our simulations, with a view to the collective nature of the excitations in ID systems, can 
be considered as constituted by both kinds of elementary excitations, appropriately weighted by the presence of the trapping 
potential. 
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We now aim at describing the dynamics that we observe numerically through the evaluation of the effective mass for the pairs 



'^doublon — 



K=0 



7- 



C/2 



16 J2' 



(33) 
(34) 



where E is defined by Eq. 



B. The frequency of the non-interacting particle in the combined harmonic trap and lattice potential 

The idea is to compare the oscillation frequency from the numerical data to the one given by the non-interacting particle in a 
combined lattice and harmonic trap potential. As a reminder, the formula for the centre of mass (COM) position was given by 



(5exp 



flho 



sin^ (y/8t) 



V 



1 



s'- . fvt 

— TT- sm — 

2«L V 4 



Neglecting the exponential prefactor time dependence, we can write the centre-of-mass oscillation frequency as 



where 



with S ~ 



(1— cos k) 



J_ 

16 



256 



n = v 

1/2 



UJCOM ~ V 



Vrricff 
, leading to 

2 1 



^COM 



lit' 



5 
aiio 



sm 



Vt 



2 



yfvfj / I - cos fc 1 
4 I mofi-y ^ 16 



VW 

256 



(35) 



(36) 



(37) 



The comparison between ujcom and the numerical data is obtained by performing a discrete Fourier transform of (n\n^ (t) for 
different values of U and N^. For = = 1, i.e. no bath, the agreement is perfect: higher interaction energies correspond 
to lower values of the oscillation frequency, in agreement with the increase of the effective mass. 
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